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A ACDM model with dark matter that decays into inert relativistic energy on a timescale longer 
than the Hubble time will produce an expansion history that can be misinterpreted as stable dark 
matter with time-varying dark energy. We calculate the corresponding spurious equation of state 
parameter, w^, as a function of redshift, and show that the evolution of depends strongly on the 
assumed value of the dark matter density, erroneously taken to scale as a"^. Depending on the latter, 
one can obtain models that mimic quintessence {w^ > — 1), phantom models (w^ < — 1) or models 
in which the equation of state parameter crosses the phantom divide, evolving from > — 1 at 
high redshift to < — 1 at low redshift. All of these models generically converge toward w,f, « — 1 
at the present. The degeneracy between the ACDM model with decaying dark matter and the 
corresponding spurious quintessence model is broken by the growth of density perturbations. 



I. INTRODUCTION 

Cosmological data from a wide range of sources in- 
cluding type la supernovae []|, the cosmic microwave 
background m, baryon acoustic oscillations |3|, cluster 
gas fractions [3, [1] and gamma ray bursts [1, Q seem to 
indicate that at least 70% of the energy density in the 
universe is in the form of an exotic, negative-pressure 
component, called dark energy. (See Ref. [Sj for a com- 
prehensive review). A parameter of considerable impor- 
tance is the equation of state (EoS) of the dark energy 
component, defined as the ratio of its pressure to its den- 
sity: 



DE- 



(1) 



If the dark energy is due to a cosmological constant, 
then w is constant and exactly equal to —1. A value 
of w = — 1 is consistent with current observations (9l.[l0|. 
On the other hand, a variety of models have been pro- 
posed in which w is time varying. Perhaps the simplest 
model for time- varying w is to take the dark ener gy t o be 
due to a scalar field, dubbed "quintessence" |lll - [l9| . As 
the data continue to improve, it is important to be able 
to distinguish models (such as quintessence) with a true 
time- varying w from ACDM models which can effectively 
mimic such a time variation. 

To motivate our investigation, consider first the simple 
ACDM model, with present-day dark matter and cos- 
mological constant densities of pdmo and p\o, respec- 
tively. (Zero subscripts will be used throughout to de- 
note present-day values of cosmological quantities). If 
one does not know a priori the value of pdmo ot the 
fact that the dark energy is a pure cosmological constant, 
then part of the dark matter density can be absorbed into 
the dark energy, producing (erroneously) a time-varying 
dark energy. In particular, the ACDM model with pdmq, 
PAo is degenerate with the quintessence model having a 
present-day dark matter density pdmo Pdmo, and a 
time-varying dark energy component with density p^ = 
PA + Pdm — Pdai and equation of state 



PAO 



Pao + {pdmo - Pdmo)(1 + zY- 



(2) 



(Since our paper deals with spurious measurements of a 
time- varying equation of state, we will use tildes through- 
out to deal with spurious/unphysical/mismeasured quan- 
tities; non-tilded variables will refer to true physical 
quantities). This degeneracy has been exhaustively ex- 
plored in Refs. |20l423| . who pointed out that it cannot 
be resolved without independent knowledge of the dark 
matter density. 

On the other hand, if observational data led to an 
equation of state for the dark energy that mimicked the 
evolution of given in Eq. ([2]), Occam's razor would 
suggest that the correct model for the universe was ac- 
tually ACDM, with the appropriately different value for 
the dark matter density. It is of interest to determine 
if there are any less obvious spurious time- varying equa- 
tions of state for the dark energy that arise from simple 
variations on ACDM. 

One such model is ACDM with decaying dark matter. 
If the dark matter decays into inert radiation, and one 
interprets the measured behavior of H{z) under the in- 
correct assumption that the dark matter is stable, then 
the best fit to the observations will be a quintessence 
model with a time- varying equation of state. This occurs 
because the nonstandard time- variation in both the dark 
matter and inert radiation densities gets absorbed into 
time variation of the dark energy. This was first noted 
by Ziaeepour [l^], who argued that it leads, in general, 
to a mismeasured value of the equation of state param- 
eter satisfying < —1. (Note that similar < — 1 
behavior occurs when energy can be exchanged between 
dark matter and dark energy (2^, but this is a different 
class of models than those discussed here). 

In this paper we reexamine the behavior of such mod- 
els, and provide an improved calculation of w^{z). We 
show that the behavior of w^{z) is very sensitive to the 
assumed value of pdm- Different choices for pdm can 
lead to models that mimic quintessence, phantom-like 
models, or models that cross the phantom divide. In the 
next section, we present the calculation of w^{z), and we 
discuss our results in Sec. III. 
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II. THE EFFECTIVE EQUATION OF STATE 

We consider a ACDM Universe which consists of bary- 
onic matter with density ps, dark matter with density 
Pdmi and cosmological constant with density pA- The 
dark matter is assumed to decay at a rate T into invis- 
ible relativistic energy with density p^. By "invisible", 
we mean that the relativistic decay products are assumed 
to interact very weakly with ordinary matter. If the de- 
cay products did interact, e.g., electromagnetically, they 
would be easily detectable at the lifetimes of interest here, 
and the model would already be ruled out. 

The equations for the evolution of matter and the 
decay-produced radiation are the following: 

RDM = —'iHpDM — TpDM, (3) 

PH = -AHpji + TpDM, (4) 
Pb = -3HpB, (5) 

where H is the Hubble parameter given by the Friedman 
equation 

3i?^ = 8ttG {pdm + Pb + Pr + Pa) , (6) 

and dots denote time derivatives. We ignore "ordinary" 
radiation, since it has a negligible effect on the expansion 
rate in the redshift regime relevant to this problem. 

The existence of such decays with lifetimes comparable 
to the Hubble time can be constrained by their effects on 
the CMB, large-scale structure, and SN la observations 
[26l - l29j . We will confine our attention to lifetimes satis- 
fying the constraint given in Ref. [1^: F^^ > 100 Gyr, 
or rto < 0.15, where to is the present-day age of the 
universe. 

The PR term in the expression for H, as well as the 
nonstandard evolution of rdm, will clearly lead to an 
effective equation of state different from the ACDM cos- 
mology. Observers unaware of the decaying nature of the 
dark matter might try to explain the time- varying equa- 
tion of state with the help of a quintessence field (f) in 
addition to ordinary non-decaying dark matter. In their 
model, (denoting the density of presumed non-decaying 
dark matter by rdm) the same expansion rate H will be 
given by: 



3i?^ = 8nG {pDM + Rb + Ptj,) 



(7) 



Equating Eq. Q and Eq. ([T]), one can readily deduce 
the energy density of this fictitious scalar field in terms 
of the densities of matter and decay-produced radiation 
as: 



= PA + PR + RDM - RDM- 



(8) 



From this, it is straightforward to determine the effective 
equation of state of this fictitious dark energy component 
from the relation —3 (1 + w^) = din p^/ din a (and also 
using equations ^ and (g])): 



W0 



Pfl/3 - PA 



PA + PR + RDM - RDM 



(9) 



Note that does not depend explicitly on ps- However, 
there is an implicit dependence, since we assume a fiat 
geometry with f2 = 1 , so that floM + + ^r. = 1 — ^b- 
(In this paper, we will use ^dm, ^^a,--- to refer only to 
present-day values, so we drop the zero subscript in these 
cases) . 

Consider first the qualitative behavior of w^. The so- 
lutions to Eqs. ([3]) and (|3]) can be written in terms of the 
scale factor a(t) as [lO] 



RDM = Rdmq I — I exp [-F(i - to)] (10) 



PR = Rdmo I — ) X 



V 



exp (-Ft') rdt', (11) 



while the fictitious non-decaying dark matter density 
evolves, of course, as 



(12) 



Pdm = Rdmo — 



Clearly, the behavior of depends on the assumed 
value of Rdmo- Using SN la data alone (or more gen- 
erally, the behavior of H{z) alone) there is no best-fit 
value for the dark matter density, once the equation of 
state of the dark energy is allowed to be a free func- 
tion of z, just as in the case of non-decaying dark matter 
discussed earlier [20l - [23| . and the best one can do is to 
derive the behavior of u;^ as a function of Pdmq- Of 
course, there are other cosmological measurements that 
constrain the dark matter density. For instance, large- 
scale structure provides a constraint on the present-day 
dark matter density, while the CMB constrains the dark 
matter density at high redshift, but taken together these 
limits are consistent with a change in the comoving dark 
matter density as large as 15%, as noted previously [29j . 
However, it is reasonable to assume that an observer er- 
roneously postulating stable dark matter would derive a 
value for rdmo that lies somewhere between the values 
obtained by taking rdm = Rdm at high redshift (Ft ^ 1) 
or by taking rdm = Rdm at the present day (t — to). 
Equating rdm from Eq. pO)) and rdm from Eq. (fT2)) 
at t — > and t = to then gives the plausible bounds on 
Rdmo- 



Rdmo < Rdmo < Pdmo exp(Fto), 



(13) 



where the upper bound corresponds to equality between 
RDM and RDM at high redshift, while the lower bound 
assumes equality today. 

To parametrize this uncertainty, we will take 



Rdmo = Rdmo exp(AFto), 



(14) 



where < A < 1. Here A = corresponds to setting 
the (spurious) stable dark matter density equal to the 
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decaying dark matter density at the present, while A = 1 
corresponds to equaUty t — >■ 0. 

To gain some quahtative insight into the behavior of 
w^, it is necessary to derive an approximation to the 
decay-produced radiation density in Eq. ([TT|) . Note that 
we always have exp(— Ft) « 1 — Ft for our observational 
bound of F^o < 0.15. This approximation can be used to 
integrate Eq. pT|) in the matter dominated era (z ^ 1) 
to yield [13 



PR = -TtpDMoia/ao) 
5 



(15) 



This result becomes progressively less accurate as the 
cosmological constant begins to dominate at late times, 
but it will be sufficient to provide qualitative insight into 
the evolution of w^. 

Substituting the expressions for pdm, Pdm and pu 
from Eqs. and (HH), respectively, into Eq. ©, 

and expanding out to lowest order in Tt and TIq (both 
of which are ^ 1), we obtain 

~ -pA + lrtpDMo{l + zf 

PA + [(i-A)rto-|rt]pi,Mo(i + ^)3' ^ ' 

In order to express as a function of redshift alone, 
we make the approximation that the expansion law never 
diverges very far from ACDM (see, e.g., Ref. (H). In 
this case, the cosmic time t is related to the redshift as 



t{z) = sinh 



■[r!A/(l-l]A)]l/2 



(1 + ^) 



3/2 



(17) 



where t\ is a constant with units of time that is related 
to the energy density of the cosmological constant by 



tA = V1/(6^Gpa). 



(18) 



Eq. PT| will be a good approximation as long as 
is close to —1 whenever the dark energy dominates the 
expansion. In the models we investigate here, can 
significantly diverge from —1 only at early times, when 
the expansion is dark-matter dominated, so we expect 
Eq. (|17p to be sufficiently accurate for our purposes. 
Then Eq. can be written as 



-»A + (l/5)F^o/(z)(l - nA){l + zf 
flA + [(1 - A) - (2/5)/(z)]Fto(l - r!A)(l + zf ■ 



(19) 



where f{z) is the fractional age of the universe at redshift 
z, given by 



t{z) 
to 



sinh ^ 



[r!A/(l-f^A)]'/'(l + 2)-3/'' 



tanh 



. (20) 

In Figs. 1-2, we show the exact evolution of w^, de- 
rived by numerical integration of Eqs. (jS]) and (U), along 



with our approximation given by Eqs. (|T9)) and ([20]) . for 
different values of Tto, A, and ft a- We see that our ana- 
lytic approximation is quite accurate for all of the cases 
examined. 



III. DISCUSSION 

Eqs. (|19p and ([20)1 yield several important insights into 
the behavior of w^. First, we note that all of our models 
converge to — > — 1 at late times. This is obvious from 
the fact that pA is always the dominant component at 
late times. However, the way in which evolves with 
redshift is extremely sensitive to the assumed value of 
Pdm through its dependence on A. For small values of 
A (i.e., taking pdm nearly equal to the present-day value 
of Pdm) the spurious time- varying dark energy evolves 
as a quintessence component with > —1, with 
decreasing with time. In the terminology of Ref. [sH. 
the model mimics a "freezing" quintessence field. 

For larger values of A, the model crosses the phantom 
divide at some redshift Zc, with u;^ > —1 at z > Zc and 
< —1 at z < Zc- An approximate value for Zc can be 
derived from Eq. 



1 



Qa 



(1 - 17a) sinh^[5(l - A) tanh" 



1/3 



(21) 



Note that Zc is independent of F (which is apparent in 
Figs. 1-2) and depends only on A. The expression for 
the time tc at which crosses —1 can be derived from 
Eq. (1161) and takes the particularly simple form 



ieAo = 5(l-A). 



(22) 



From Eq. ([^^ . it is obvious that these models cross the 
phantom divide whenever 0.8 < A < 1. (Strictly speak- 
ing, models with A < 0.8 also cross the phantom divide, 
but this crossing takes place in the future, while A = 0.8 
gives exactly equal to —1 at the present). Finally, for 
A = 1 (i.e., setting pdm equal to pdm at early times) 
the model behaves like a phantom, with < — 1 at all 
times. 

The ACDM model with decaying dark matter can be 
distinguished from the corresponding quintessence model 
with the same H{z) using the growth of density pertur- 
bations. (For a general discussion of linear perturbation 
growth as a probe of dark energy, see, e.g., Refs. [Fi,[32|- 
|42||). The equation for the evolution of the density per- 
turbation, S, in the linear regime {S <^ 1), well inside the 
horizon, is 



S + 2Hd- AttG{pdm + Pb)S - 0, 



(23) 



where the dot denotes the derivative with respect to time, 
and we assume the baryons are decoupled from the pho- 
ton background and can cluster freely. By construction, 
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FIG. 1; The evolution of the spurious equation of state parameter, w^, in a ACDM model (JIa = 0.7) with decaying dark 
matter, in which the dark matter is erroneously taken to be stable and the corresponding time variation is absorbed into a 
time-varying equation of state for the dark energy. Here V is the decay rate of the dark matter, to is the age of the universe, 
and A is determined by the assumed value of the dark matter density, erroneously taken to scale as d~'^ . Broken lines denote 
the exact evolution and solid lines denote the analytic approximation given by Eqs. (|19l) - (|20[) . 



the value of H{z) is the same in the ACDM model with 
decaying dark matter and in the corresponding spuri- 
ous quintessence model with stable dark matter. How- 
ever, these two models differ in the value of pdm, which 
drives the growth of density perturbations in equation 
(1^51) . In the ACDM model with decaying dark matter, 
Pdm is given by Eq. pU)). while the corresponding spu- 



rious quintessence model has a (spurious) matter density 
PDA4 given by Eq. p^ . Then using equation we 
see that 

= exp[rt + ( A - l)Ro] . (24) 
The growth rate in the spurious quintessence model will 
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FIG. 2: As Fig. 1, with ^a = 0.75. 



be larger (smaller) than the corresponding growth rate in 
the ACDM model as Tt-I- (A - l)rto is greater than (less 
than) zero. In the limiting case A = 1 the growth rate 
will be larger for the spurious quintessence model than for 
the decaying dark matter model over the entire evolution 
history, while the reverse is true for A = 0. Independent 
of A, the ratio of the perturbation growth rate at late 
times to the growth rate at early times is smaller for the 
decaying dark matter model than for the corresponding 
spurious quintessence model with the same H{z). This 



result follows from the fact that the dark matter density 
decreases more rapidly with scale factor in the former 
model. 

An obvious question is whether or not these results 
can ever be relevant. If dark matter really were unstable, 
with a lifetime much longer than the age of the universe, 
then it is certainly plausible that a cosmological signa- 
ture of unstable dark matter would be detected first in 
some combination of CMB and large-scale structure ob- 
servations. In this case, the information on decaying dark 



6 



matter would simply be incorporated into calculations for 
the equation of state of the dark energy, rendering our 
calculations moot. On the other hand, it is also plausi- 
ble that such effects would be detected first in precision 
measurements of the dark energy equation of state, in 
which case our expressions for provide a useful guide 
to the sort of spurious signal produced by decaying dark 
matter. 

These results can be generalized, in an obvious way, to 
quintessence models (or other dark energy models with a 
time-varying equation of state) with unstable dark mat- 



ter, to derive the corresponding change in the measured 
equation of state parameter. Such models, however, must 
be considered somewhat baroque. 
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